Wave function Monte Carlo method for polariton condensates 
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We present a quantum jump approach to describe coupled quantum and classical systems in 
the context of Bose-Einstein condensation in the solid state. In our formalism, the excitonic gain 
medium is described by classical rate equations, while the polariton modes are described fully 
quantum mechanically. We show the equivalence of our method with a master equation approach. 
As an application, we compute the linewidth of a single mode polariton condensate. Both the 
line broadening due to the interactions between polaritons and the interactions with the reservoir 
excitons is taken into account. 

PACS numbers: 03.75.Kk, 67.90.+Z, 71.36. +c. 
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I. INTRODUCTION 



The achievement of strong light-matter coupling in 
semiconductor microi and nanocavities^i^ has led to the 
manifestation of several novel quantum effects in the solid 
state. A prominent example is Bose-Einstein condensa- 
tion of exciton-polaritons^. A second important topic 
is the dynamics of nanocavities in the strong coupling 
regime with a quantum dot, where lasing was recently ob- 
served^. These phenomena are characterized by a large 
degree of complexity, due to the coupling of quantum 
dynamics with the nontrivial dynamics of the injected 
incoherent particles. The physics that these solid state 
systems have in common is that there are some degrees of 
freedom for which the quantum mechanical nature is im- 
portant, while other degrees of freedom are essentially 
classical. This is illustrated in Fig. [1] Due to their 
coupling, it is not straightforward to construct a model 
for the combined dynamics of the classical and quantum 
parts. The typical approach is an adiabatic elimination 
of the classical reservoir, assuming that it adapts instan- 
taneously to the state of the quantum parli^. For polari- 
ton condensation, this type of approach was developed 
by Laussy et al.— . Also in Keldysh Green function treat- 
ments, the bath degrees of freedom are usually integrated 
out^, which again limits the types of interactions between 
bath and system that can be taken into account. 

In this Article, we want to construct a theory that 
goes beyond such approximations and does not make any 
other approximations apart from a classical treatment of 
the reservoir and that the coupling between the reser- 
voir and quantum system is weak. We will show that 
a master equation approach leads in general to an in- 
finite hierachy of equations, that does not allow for an 
easy truncation. We subsequently show that a quan- 
tum jump model can be constructed that is equivalent 
to the infinite set of coupled master equations. We ex- 
plicitly show that the Schawlow-Townes linewidth and 
Henry linewidth enhancement factor^'^*^ and line broad- 
ening due to polariton-polariton intcractions^i are repro- 
duced with our model. 




FIG. 1: Sketch of the problem under study: The incoherent 
reservoir particles relax at a rate r^ to the quantum sytem by 
giving the excess energy to the environment modes p. 



II. MASTER EQUATION 

We will derive our equations for the quantum evolution 
of the coupled polariton-reservoir system, assuming the 
interaction between the two to be of the form 



Hi = ^ gkqp\bla + h.c. 



(1) 



This Hamiltonian describes the relaxation process illus- 
trated in Fig. [TJ a reservoir particle bk relaxes to the 
quantum system mode a, depositing its excess energy in 
the 'phonon' mode Pg, that we assume for simplicity to be 
at zero temperature. In the usual Born-Markov approxi- 
mation, the following master equation can be derived for 
the reduced density matrix 



9 



dt 



—Tt:r{p) = ^ -^Tr_R, ia^bkpb\a - b\bkaa} p + h.i 



.(2) 
under the assumption that coherence in the reservoir is 
neghgible, so that e.g. Tr/j[6^(a''')^p] = 0. The rate r^ 
is the Fermi golden rule transition rate from the state bk 
into states a and Pq . 

Unless the reservoir density is assumed to be not af- 
fected by the polariton dynamics, Eq. ([2]) is not a closed 
equation of motion for the polariton density matrix ps = 
Trfl (p) , but it is coupled to the reservoir weighted density 
matrices pk = Trji{b[,bkp) as 



d_ 
dt 



k 



(3) 



The dynamics of pk is under the same assumptions 
that lead to Eq. ^ described by 



d_ 
dt 



Pk 



Eli 
2 
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{aa) pj,k + P],kaa^) + Vj 



(4) 



where we have introduced pk.j = TT:B,{h\,hkh-bjp). 

For the pumping term V in Eq. (j?]), we take a stan- 
dard gain mechanism that does not include gain satura- 
tion (the physical pump is assumed to be for example an 
external laser that excites the reservoir), leading to 

[b\pbj+bjh]p + \\.c. 



V, 



Tr 



b\b. 



j 3 2 
PTrR (p) . 



(5) 
(6) 



For the last step, we have assumed that the occupation 
of each mode is on average much lower than one, so that 
we can neglect Ti{bJ^bkp) whith respect to Tr(p). 



III. QUANTUM JUMP MODEL 

It is clear that the dynamics of pj,k will involve the ex- 
pectation value of densities in three reservoir modes and 
so on. Physically, Eq. ([2|) expresses that gain is propor- 
tional to the density matrix weighted with the reservoir 
densities. For weak gain saturation, i.e. when the num- 
ber of reservoir particles weakly depends on the polari- 
ton number, the gain saturation can be treated pertur- 
batively^. It is the aim of this article to go beyond such 
a perturbative treatment and propose a wave function 
Monte Carlo method (see e.g.J^) that is equivalent to the 
hierarchy of Eqns. ©, (g]), etc. 

Because the reservoir is assumed to be classical, its 
state can be described by the occupations alone r — 
{ni, 712, • ■ •}• The polariton part is treated fully quan- 
tum mechanically by a wave function 1$). The state S of 
the system and reservoir is then S = [r, |$)]. A stochastic 
dynamics for S leads to a probability distribution P{S), 
that has an associated density matrix 



p^ d<^ 



ri,r2,... 



P({nJ,|$))|$)($|, 



(7) 



where the integration runs over all possible values of the 
wave function. The reservoir weighted density matrices 
are defined analogously as 



Pk = / d^ 



ri,r2,... 



P({n,;},|$))nfe|<i>)($|. 



(8) 



For the deterministic part of the dynamics, we propose 
the nonlinear dynamics that is standard in wave function 
Monte Carlo: the reservoir densities are invariant under 
the deterministic part of the evolution drii/dt — 0, where 
the wave function evolves according to 

■^ \^) ^^Y.^knkaa^\>^) + \R{{n,})\\a^\m^W), 

k 

(9) 



'dV 



where R{{ni}) = ^^ r^nk is the total spontaneous rate 
for the creation of polaritons out of the reservoir. The 
second term in Eq. (|9]) compensates for the norm- 
reducing dynamics due to the first term. After an evolu- 
tion over time 5t, the wave function is transformed to 



mt + 5t)) 



1 



VT^^ 



\m)-^^^aa^\^) 



where 



5p^5tR{ni)\\a)\<^)f 



(10) 



(11) 



is the probability that a polariton is created out of the 
reservoir during the time 5t. 

The jump part of the dynamics due to the creation of 
a polariton from a reservoir particle in mode k is given 

by 



Uk ^ Uk- 1, 
at|$) 



m 



|aT|$) 



(12) 
(13) 



Such a jump occurs in an interval dt with probabilities 
5pk = 5trknk\\a'\^)\\'^ that sum as J^k^Pk = ^P- 

Under the dynamics (|T0)) . (|12l) and (TTSl) . the expecta- 
tion value of the density matrix |$)($| evolves as 

E[mt + dt))mt + st)\] = E[mt)){<^{t)\ 



--(1- Sp)R(nMt^^}^^ 



1 



(1 - 5p)R[n^)5t 



1$) 



5p 



op 



— Sp 
)($|a 



^/WsiVWSi 



(14) 



With the definitions (g]) and ([TTl), Eq. dH]) implies that 
the expectation value E[\^{t + St))(^{t + St)\] obeys the 
equation of motion ^. 

In the same way, one can verify that the evolution 
stochastic evolution for | $) leads to the equation of mo- 
tion (|31) with the definition 



p,,k^Jd<^ E P{{n^},m)n,nkm{^\ 



(15) 



The correspondence between the master equation de- 
scription and the quantum jump appoach is entirely clear 
from a physical point of view. The classical reservoir con- 
tinuously monitors the polariton mode, which leads to 
the quantum jumps. The present approach is a general- 
ization of the usual quantum jump models in the sense 
that we take into account the backaction of the system on 
the reservoir, which in turn affects the system dynamics. 

The pumping of the reservoir can be taken into account 
by the additional jump process acting on the reservoir 
state alone Uk — > n^ -I- 1 with a probability in the interval 



6t equal to 6p = 6tP. This jump process leads to a term 
in the evolution of the density matrix p^. of the form 



^pfc = ^/rf$E^^i*)('^i' 



(16) 



which reproduces Eq. ([6]). In addition, all internal dy- 
namics in the reservoir, due to for example collisions be- 
tween particles in the reservoir can be simply modeled 
by applying the modification on the reservoir mode oc- 
cupations and keeping the polariton state the same. 

The dynamics ((9|l- (fT3|) naturally leads to polariton 
states that arc in a number state, because the reservoir 
monitors the polariton number by the stimulated relax- 
ation. This implies that a semiclassical stochastic model 
for a quantum dot embedded in a nanocavity, such as 
e.g. performed in Ref J^ does not make any approxima- 
tions apart from the ones involved in the derivation of the 
quantum Markov dynamics ^. Moreover, this picture 
is consistent with the discussion of the coherence proper- 
ties of laser light by M0lmer'^^, who advocates the point 
of view that the photon state inside a laser cavity is a 
number state, rather than a coherent state. 



IV. TEMPORAL COHERENCE 

For the computation of the temporal coherences, the 
stochastic wave function dynamics can be computed on 
a doubled Hilbert spaced for the polariton state. More 
specifically, the first order coherence g^ \t -\- T^t) = 
{a){t + T)a{t)) can be computed by propagating the state 
Se = [r, 0], where the normalized state vector 9{t) is pro- 
portional to [$(i) a$(i)]"^ over a time r according to the 
prescription detailed in Ref J^, that is to apply the evo- 
lutions ^ and (TTSl) on the two components of the wave 
function with a proper normalization of the total wave 
vector. The temporal coherence is then given by the ex- 
pectation value 



.g(i)(t+r,i) =i?[|| {'^{t)\a\'^{t))\\\e^{t+sr)\c 



t| 



-ST)\a'\t)l{t + T))] 

(17) 
, where 9i and 62 refer to the first and second part of the 
doubled Hilbert space respectively. 

To understand the coherence properties of the polari- 
ton field, it is instructive to compute the temporal coher- 
ence of a state that is initially in a number state with a 
large number of polaritons N that decay at a rate 7. The 
deterministic evolution of the state = {u\N),v\N —1))'^ 
gives a change of the components on the doubled Hilbert 
space 



d_ ( u\N) 
dt UlA^-l) 






(18) 



The coherence is proportional to uv = u\/l — u^. When 
a polariton is lost due to a quantum jump after a time r. 



the change in u under a deterministic evolution followed 
by quantum jump is 



u -)■ 



yjnv? + [N — 1)(1 — V?) exp[rT] 



(19) 



Under this evolution, u tends to zero for long times. The 
temporal coherence time is thus equal to the decay time 
of u for u tending to zero. For small u, the waiting time 
distribution of r is given by 



P(t) = (A^-l)7exp[-7(7V- l)r]. 



(20) 



Averaging the evolution of u over this waiting time in 
the limit for small u and to leading order in 1/A^ leads to 
M — >■ [1 — \/{%N'^)]u. The average number of jumps per 
unit time is ^N , so that on average 



du 

It 



8N 



(21) 



The same calculation can be repeated including the 
gain due to the reservoir r. The change in u for a sin- 
gle quantum jump, either loss or gain, is still given by 
Eq. p^ in the large N limit. The number of jumps how- 
ever doubles (on average the gain r should compensate 
for the losses 7), so that we obtain the Schawlow-Townes 
expression for the decay of the coherence: 



7c 



7 

'an' 



(22) 



The equal contribution of losses and gain to the laser 
linewidth was noted by Scully and Lanib^'^. 

The equation (j22p is entirely general and does not make 
any assumptions on the nature of the gain medium. In 
the present formulation, the robustness of the Schawlow- 
Townes linewidth comes from the fact that the dccoher- 
ence per quantum jump does not depend on 7 or r. It is 
therefore insensitive to the fluctuations in the gain and 
the strength of gain saturation. 

The physics becomes more complicated when inter- 
actions between the reservoir and polaritons are taken 
into account. This is important in semiconductor lasers 
where the carrier concentration affects the refractive in- 
dex and thus shifts the polariton energy. Fluctuations in 
the reservoir occupancy then affect the laser line width, 
leading to the well known Henry linewidth enhancement 
factor^S. We will now proceed to discuss the enhancement 
of the linewidth due to such interactions. They result in 
an additional phase shift of u under the deterministic 
evolution, that is taken to be of the form A^ = gnAt. 

A situation in which the linewidth enhancement can 
be easily computed is the one illustrated in Fig. [51 Here, 
the gain due to a single particle in the reservoir is as- 
sumed to be much larger than the cavity line width. As 
a consequence, the reservoir occupation has a very small 
probability to be occupied by more than one particle. 
The evolution of u under a single gain peak is given by 



u -^ 



\fne 



^nv? + (N - 1)(1 - m2) exp[7T] 



(23) 



iUUi 



FIG. 2: Gain as a function of time in the case of a strong 
backaction of the cavity dynamics on the reservoir population. 



Averaging again over the waiting time distribution for r, 
and adding the decoherence due to the losses, we obtain 
in the large A^ limit 



d 



.9r 



1 + 4(to/r)2 



47V 



u, 



that gives a polariton line width 



7c = 



1 + 4(W0' 



4N 



(24) 



(25) 



The 1/A^ scaling with the number of polaritons is of 
the Shawlow-Townes form (1^^ and the increase of the 
linewidth is the Henry linewidth enhancement^. The 
standard deviation of this factor is in a semiclassical ap- 
proximation, that assumes small fluctuations in the gain 
medium density^*^ and is limited to large polariton oc- 
cupation numbers. It thus turns out that the Henry 
linewidth enhancement is also a robust phenomenon that 
is insensitive to the details of the system-reservoir cou- 
pling. 

In our formalism, this semiclassical approximation cor- 
responds to approximating the quantum jump model by 
a diffusion modcl^^. The Langevin equations for the pho- 
ton and reservoir occupations Nr are taken to be 



dN = [rNniN + 1) - -/N]dt + y^rNR{N + l)dWG 

+ y/^dWL, (26) 

dNR = -rNR{N + l)dt-^/7N^JWTT)dWG. (27) 

Here dWcL are independent Gaussian noise terms with 
variance {dWcLdWcL = dt). The stochastic terms rep- 
resent the shot noise fluctuations of the photon gain and 
losses. 

In the large photon limit, the Langevin equation (P7)) 
can be solved approximately by linearizing around the 
mean field solution {N'^^\N^^^) : N = N^°^ + 5N and 
Nr = N^°^ + 5Nr, yielding 



d 
di 



SN 
SNr 

1 
-1 



rNo 

-7 -rNo 



SN 

SNr 



y/rNR^oNodWa 



1 



^WodWL. (28) 



The interaction with the reservoir induce a phase shift 
between the two components of in Eq. (|18p , that is equal 
to 



<Pit) = CJR / SNR{t')dt'. 

Jo 



(29) 



This phase shift causes a decoherence that is to be 
added to the Shawlow-Townes decoherence mechanism 
that leads to a decrease of the absolute value of the first 
component of the vector Eq. (|18l) . 

Using the expression for Gaussian variables 

first order coherence as 



,(1) 



5:e^-poiW=exphA(i)], 



where 




t -|2\ 

SNR{t')dt' 



(30) 



(31) 



Here, the coupling constant qr quantifies the strength of 
interactions between reservoir excitons and polaritons. 
The subcript in Eq. (|30|) indicates that we only con- 
sider the effect of the rcservoir-polariton interactions on 
the polariton linewidth. Other contributions to the de- 
coherence, such as the Shawlow-Townes mechanims and 
the effect of polariton-polariton interactions (see below) 
have to be added to this decoherence. 

The expectation value in the above expression can be 
rewritten as^^ 



n2\ 



5NR{t')dt' 







= 2 / dt'{t - t')$fl(t'), (32) 







where we have used that the reservoir correlation func- 
tion 



^R{ti - h) ^ {SNR{h)5NR{t2)), 



(33) 



only depends on the time difference. 

Eq. p2|) is a relation that is used in used in the stan- 
dard Kubo theory for the line width of an emitter with 
time-dependent frequenc y^^'^^ . Two limiting cases of this 
expression are of particular interest. For times that are 
short compared to the correlation time of the reservoir 
fluctuations tr, the coherence decreases as a Gaussian 



X{t) ^ $fl(0)t2 



(34) 



where for times much larger than the reservoir fluctuation 
time 



X{t) ^ $fl(0) TRt « $fl(0) f 



(35) 



The above inequality is physically described as motional 
narrowing: for times much longer than the fluctuation 
time of the reservoir, the decoherence effect on the system 
is suppressed. 



The reservoir correlation function can be easily com- 
puted in the Fourrier domain, using the relation 



^R{t) 



duj 
2^' 



'{\SNr{u)\' 



(36) 



From the Fourrier transform of Eq. (|28| , one obtains 

LO^ + 2/72 



where the linear damping rates are given by 
ri,2 = i (rNn T \/(riVo)2 - A-/rNo. 



Whittakerli. Note that it scales very differently than 
the linewidth induced by the interactions with the reser- 
voir (p9l) . In particular, the line-broadening due to the 
polariton-reservoir interactions decreases with increas- 
ing polariton number, where the broadening due to the 
polariton-polariton interactions increases with increasing 
polariton number. 

For short times, the coherence decays as 



(37) 



(38) 



Ji) 



.9poi-poi(*) = exp 



, (27 + rJVo) 
■9p 4r 



Using Eqns. (pO| . (pGj) and ([37|). the motional nar- 
rowing dominated long time behavior of the decoherence 
induced by the interactions between reservoir particles 
and the system polaritons is given by 



Ji) 

71CS— pol 



-Doi(0 '- exp 



m 



Nn 



T 



for t>i/ri, 

(39) 

which is despite the very different approximations, iden- 
titical to the interaction contribution in Eq. (1^^ . 

The Shawlow-Townes contribution to the coherence de- 
cay (f22|l is to be added to this expression. Mathemati- 
cally, it comes from the decrease in magnitude of u{t) in 
Eq. (fT8|) . where the decay (p9|) is due to the scrambling 
of its phase. 

The Gaussian short time behavior is governed by 



5rcs--pol(*) 



g|j7(27 + WVoy^2 



exp -^1^^1-L^T^ for T « l/T, 

. (^^^ 
The Gaussian early time decay was not obtained in the 

case of the large gain saturation, described by Eq. ([24| . 
This is consistent with the fact that the reservoir corre- 
lation time vanishes in this approximation. 

Entirely analogously, the line broadening due to the 
interactions between the polaritons can be computed. 
Within the linearized model p8|) . the photon fluctuations 
are 



{\SN{co)\' 



2a;2 + (riVo)^ 



(c.2 + rf)(c.2+ri) 
The contribution to the long-time decoherence is 



(41) 



9pol-pol(^) 



exp 



-^-27^ 



for r>l/ri, (42) 



where gp is the interaction strength between polaritons. 
The dependence on particle number Nq and linewidth 
7 is the same as the one obtained by Eastham and 



for T < l/ri,2. 

(43) 

Compared to the result for reservoir-polariton interac- 
tions PO)) . a large polariton population is less favorable 
for good temporal coherence. 

It is also interesting to notice the different dependence 
of the coherence on the gain saturation parameter r. The 
coherence times due to the reservoir-polariton interac- 
tions in Eqns. ([5^ . (^(71 scale as ^ r^ . The coherence 
time in Eq. ((321) scales weaker, as ~ r, and at long times, 
the coherence decay (j39l) even becomes independent of 
the gain saturation. 



V. CONCLUSIONS 

We have constructed a quantum jump model that is 
able to describe the dynamics of a coupled quantum- 
classical system, where the quantum system enjoys gain 
from the classical one. Our model makes it possible to 
describe arbitrarily complex dynamics in the reservoir 
and energy dependent gain. As a first application of 
'2 'the theory, we have investigated the linewidth of a sin- 
gle mode polariton condensate. The Shawlow-Townes 
line width, line broadening due to polariton-reservoir in- 
teractions and polariton-polariton intereactions can be 
described with this formalism. We have shown that the 
these two mechanisms for line broadening exhibit a dif- 
ferent dependence on the number of polaritons and on 
the reservoir gain saturation. 

Our approach can be staightforwardly combined with 
a full Boltzmann dynamics of the reservoir, so to give an 
ab initio quantum description of polariton condensation 
in semiconductor micro and nanocavities. Moreover, the 
formalism can be applied to spatially extended systems 
as well, that require to take into account the energy de- 
pendence of the gain processi^. 
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